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Abstract: We examine in detail the neutrino masses and mixing patterns in an extension 
of the minimal supersymmetric standard model with three gauge-singlet neutrinos and It- 
parity violation. The Majorana masses for the gauge-singlet neutrinos as well as the usual 
//-term for the Higgs superfields are generated at the electroweak scale through the vacuum 
expectation values of the singlet sneutrinos. The resulting effective mass matrix for the 
three light neutrinos have contributions from the seesaw mechanism involving the singlet 
neutrinos as well as due to the mixing with the heavy neutralinos. This model is popularly 
known in the literature as the from v supersymmetric standard model" (/wSSM). We 
show that even with flavour diagonal neutrino Yukawa couplings, the global data on three- 
flavour neutrinos can be well accounted for in this scenario, at the tree level. We also 
analyze the mixing in the chargino and the Higgs sector and calculate the decays of the 
lightest supersymmetric particle in this model. The decay branching ratios show certain 
correlations with the neutrino mixing angles, which can be tested at the LHC. Some other 
phenomenological implications of such a model have been discussed. 
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1. Introduction 



The experimental results on neutrinos provide strong evidence in favour of non-zero neu- 
trino masses and mixing angles. Various neutrino oscillation experiments suggest that the 
mixing pattern of the three light neutrinos is bilarge, that is to say there are two large 
mixing angles and one small mixing angle. The data can be explained very well with the 
following set of mass squared differences and mixing angles Jll ||| 

7.05 x l(T 5 eV 2 < Am^ < 8.34 x l(T 5 eV 2 (3a), (1.1) 
2.07 x 10" 3 eV 2 < lAmlil < 2.75 x 10~ 3 eV 2 (3a), (1.2) 



30° < U < 37°, 35° < #23 < 54°, 13 < 13° (3a), (1.3) 



where Am 2 , = — 



In order to explain the results in eqs.(l.l)- (pq), one needs to go to a theory beyond 



the standard model (SM). A very interesting possibility to look for new physics is supersym- 
metry (SUSY) which has the ability to provide a solution of the so-called 'gauge hierarchy 
problem' connected with the mass of the Higgs boson. SUSY predicts new particles at the 
TeV scale which can be tested at the forthcoming Large Hadron Collider (LHC). Naturally 
it is tempting to see whether TeV scale SUSY is also one of the candidates which can 
explain the observed pattern of neutrino mass squared differences and mixing. There have 
been several proposals in recent times, which attempt to explain the experimental data on 
neutrinos in the context of a supersymmetric theory. Perhaps the most well studied class 
of models in this context is the one which includes R-parity violation ||, Q in the Minimal 
Supersymmetric Standard Model (MSSM). Neutrino mass squared differences and mixing 
angles have been calculated under various assumptions and it has been found that the 
neutrino data can be explained well when contributions to the neutrino mass matrix at the 
tree and one- loop level are considered ||]. 

On the other hand, though SUSY provides some elegant solutions to accommodate the 
experimental data on neutrinos, it has been plagued by a few urgent questions which do not 
have very satisfactory answers yet. One of them is the so-called "/z-problem" [JjJ related to 
the bilinear term fiH\H2 in the MSSM superpotential. The electroweak symmetry breaking 
requires the value of \i to be roughly of the order of a few hundreds of GeV. This requires, 
in the absence of any fine cancellation, that fx is roughly of the order of the soft scalar 
masses and both of them should be somewhere around a TeV or a few hundreds of GeV. 
Since fi parameter respects SUSY, there is no obvious reason why it should be of the same 
order as SUSY breaking soft scalar masses. This defines the "^-problem". There have 
been several attempts to address the solution to this problem and all of them requires the 
vacuum expectation value(s) (VEVs) of some additional field(s) to generate the /i-term 
after the symmetry breaking. 

One of the solutions is the next-to-minimal supersymmetric standard model (NMSSM), 
which introduces a superfield S, singlet under the SM gauge group. The \i term is absent 
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from the superpotential and it arises when the scalar component of S acquires a VEV. 
This VEV is determined in terms of the soft supersymmetry breaking terms through the 
minimization condition. If the SUSY breaking scale is close to the electroweak (EW) scale 
then the effective /i term is also of the order of the EW scale. However, as in the case 
of MSSM, the NMSSM also cannot explain the observed pattern of neutrino masses and 
mixing. 

It is also important to note in this context, that another very well known mechanism 
of generating small neutrino masses and bilarge mixing angles in a SUSY model, compat- 
ible with the experimental data, is the seesaw mechanism which introduces gauge singlet 
neutrino superfields. In such cases the MSSM superpotential contains additional terms 
involving the Yukawa couplings for neutrinos as well as Major ana masses for the gauge sin- 
glet neutrinos. In the conventional scenario, the neutrino Yukawa couplings are assumed 
to be ~ 0(1), whereas the Majorana masses for the gauge singlet neutrinos are taken 
somewhere around 10 14 GeV or so. In this case light neutrino masses as small as 10 -2 eV 
can be generated. One viable alternative to the usual seesaw mechanism is to consider the 
TeV-scale seesaw. This possibility is very interesting since it may provide a direct way to 
probe the gauge singlet neutrino mass scale at the LHC and does not need to introduce a 
very high energy scale in the theory. However, in order to generate small active neutrino 
masses one needs to consider neutrino Yukawa couplings to be of order 10~ 6 . This choice 
is reasonable since we know that the electron Yukawa coupling should also be of the order 
of 1(T 6 . 

In this work we study in details, a model of neutrino masses and mixing which intro- 
duces the gauge singlet neutrino superfields (vf) to solve the /j, problem in a way similar 
to that of NMSSM. The terms in the superpotential involving the uf include the neutrino 
Yukawa couplings, the trilinear interaction terms among the singlet neutrino superfields as 
well as a term which couples the Higgs superfields to the uf. In addition, there are corre- 
sponding soft SUSY breaking terms in the scalar potential. When the scalar components 
of get VEVs through the minimization conditions of the scalar potential, an effective \x 
term with an EW scale magnitude is generated j^j. In addition, small Majorana masses of 
the active neutrinos are generated due to the mixing with the neutralinos as well as due to 
the seesaw mechanism involving the gauge singlet neutrinos. In such a scenario, we aim to 
explain the experimental data on neutrinos discussed in the beginning. In particular, we 
show that it is possible to provide a theory of neutrino masses and mixing explaining the 
experimental data, even with a flavour diagonal neutrino Yukawa coupling matrix, without 
resort to an arbitrary flavour structure in the neutrino sector. This essentially happens be- 
cause of the mixing involved in the neutralino-neutrino (both doublet and singlet) system 
mentioned above. We perform a detailed analytical and numerical work and show that 
the three flavour neutrino data can be accommodated in such a scenario. In addition, we 
observe that in this model different neutrino mass hierarchies can also be obtained with 
correct mixing pattern, at the tree level. 

In this model the neutral Higgs bosons mix with the sneutrinos and charged Higgs 
bosons mix with the charged sleptons. The corresponding scalar, pseudoscalar and charged 
scalar mass squared matrices are enlarged and we take into account the constraints on the 
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parameters coming from the positivity of the squared scalar masses. In the fermionic 
sector, in addition to the neutralino-neutrino mixing there is mixing also between the 
charginos and the charged leptons. This can also have implications for phenomenological 
studies, particularly in the context of future colliders. Because of the mixing between the 
neutralinos and the neutrinos, the lightest neutralino, which is the lightest superparticle 
(LSP) for most of the parameter space, can have novel decay modes which can be correlated 
with the neutrino mixing pattern. This can provide some unique signatures of such a 
scenario which can be tested at the LHC. 

As we have mentioned earlier, that in order to get the correct light neutrino mass 
scales, the neutrino Yukawa couplings should be of the order of 10 -6 or so. This is because 
the TeV scale VEVs of the singlet neutrinos induce Majorana mass terms of themselves 
also at the TeV scale. Similarly, the neutralino-neutrino mixing provides correct light 
neutrino mass scales only when the Yukawa couplings of the neutrinos are somewhere close 
to that of the electron. This model has been named as the "/wSSM" in ref.Q. Thus in 
this model, an interesting proposal has been given, where the generation of small neutrino 
masses and the solution to the \i problem can be accommodated with the same set of 
gauge singlet neutrino superfields without introducing an extra singlet such as in the case 
of NMSSM. The spectrum and parameter space of this model, with three gauge singlet 
neutrino superfields, were discussed in ||. However, a detailed discussion of the issue of 
neutrino masses and bilarge neutrino mixing, in order to accommodate the three flavour 
neutrino data, has been lacking and that is what we want to provide in the present work. 
We would like to emphasize here that our analysis shows that even with flavour diagonal 
neutrino Yukawa couplings, the resulting structure of the effective mass matrix of the light 
neutrinos can explain the bilarge pattern of mixing. In addition, we explore the correlation 
between neutrino mixing and the decay pattern of the LSP in this model and discuss some 
other interesting phenomenological implications. 

Other models which address the neutrino experimental data and the /i problem are 
essentially extensions of NMSSM. One of these proposals |l0[ include both the gauge singlet 
neutrino superfields ipf) and the extra singlet (S) of the NMSSM. In this case the Majorana 
masses of the singlet neutrinos are also generated at the EW scale through the VEV of 
the scalar component of S. R-parity may be broken spontaneously and the light neutrino 
masses are generated through the seesaw at the EW scale. Another possibility is discussed 
in ref.[|l]], where the effective bilinear R-parity violating terms are generated through the 
VEV of the singlet scalar S. Naturally, in this case only one neutrino mass is generated 
at the tree level whereas the other two masses are generated at the loop level. In [12], 
two neutrino masses are generated at the tree level by including explicit bilinear R-parity 
violating terms in addition to the R-parity breaking term involving S. 

The plan of the paper is the following. We start with the description of the model in 
Sec. II by writing down the superpotential and the soft supersymmetry breaking interac- 
tion terms. We also derive the minimization equations of the scalar potential and discuss 
some issues related to the vacuum expectation values of the left sneutrinos. In Sec. Ill we 
continue our discussion of the scalar sector with a more detailed look. Sec. IV describes 
the fermionic sector of the model where neutralinos mix with both the doublet and singlet 
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neutrinos and the charginos mix with the charged leptons. In Sec. V we provide a de- 
tailed discussion of the effective mass matrix of the three light neutrinos, arising because of 
the neutralino-neutrino mixing. Analytical expressions of the mass eigenvalues and eigen- 
vectors are derived under certain conditions, using degenerate perturbation theory. We 
construct the neutrino mixing matrix and show that it is possible to have two large and 
one small mixing angles in general. A detailed numerical analysis has been performed and 
we compare our results with that obtained using the approximate analytical formulae. We 
show that for realistic parameter choices, it is possible to fit the three flavour global neu- 
trino data in this scenario, even with a flavour diagonal neutrino Yukawa coupling matrix. 
The decays of the lightest super symmetric particle are discussed in Sec. VI and it has 
been observed that certain decay branching ratios are correlated with the neutrino mixing 
angles. We make concluding remarks in Sec. VII with possible future directions of our 
work. The details of various scalar mass squared matrices and the Feynman rules for the 
LSP decay calculations are relegated to the appendices. 



2. The model and its minima 
2.1 Superpotential 

In this section we introduce the model along the lines of ref. || , discuss its basic features and 
set our notations. We introduce three gauge singlet neutrino superfields, Of (i = e, /i, r), in 
addition to the fields of the minimal supersymmetric standard model. The superpotential 
of the model is written as 

W = e ab (Y^H b 2 Qtu c , + Y l d 3 mQ\d) + Y? i H?L b i e c j + Y$ H b 2 L?ty 

- e^ufHtHl + \& k i>Zi>p%, (2.1) 

where H\ and H2 are the down-type and up-type Higgs superfields, respectively. The Qi 
are doublet quark superfields, 11? [dj] are singlet up-type [down-type] quark superfields. The 
Li are the doublet lepton superfields, and the ij are the singlet charged lepton superfields. 
Here a, b are SU(2) indices, and €12 = -£21 = 1- Note that the usual bilinear fi-term 
of the MSSM is absent from the superpotential whereas two additional trilinear terms 
are introduced involving the Higgs superfields, Hi and H2, and the three gauge singlet 
neutrino superfields, vf. This is done by imposing a Z3 symmetry which is also used in 
the case of NMSSM. If the scalar potential of the model is such that non-zero vacuum 
expectation values of the scalar components (Df) of the singlet superfields i>f are induced, 
an effective bilinear term [iHfH^ is generated, where the coefficient /x = A'(P?). In the 
presence of soft supersymmetry breaking, it is usually expected that the VEVs of uf are at 
the electroweak scale. Hence the value of fi is of the order of the electroweak scale as long 
as the dimensionless couplings Aj are ~ 0(1). This gives us a solution to the so-called "/i- 
problem" . The last term in the superpotential with the coefficient k 1 ^ is included in order 



to avoid an unacceptable axion associated to the breaking of a global U(l) symmetry [13]. 
This term generates effective Majorana masses for the singlet neutrinos at the electroweak 
scale. 
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The last two terms in (|2.l| ) explicitly break lepton number (L) and hence the R-parity, 
which is defined by R = (_i )^+ 3B + 2s . Here B is the baryon number and s is the spin. Note 
that R = +1 for particles and —1 for superpartners. Since lepton number is explicitely 
broken, no unwanted massless Majoron appears in this model. One should also notice 
that the term in the superpotential involving the neutrino Yukawa couplings Y„ J , generate 
effective bilinear R-parity violating interactions e l H2Li. Here e l is determined in terms of 
the VEVs of uf and is given by e l = Yy J R-parity breaking implies that the lightest 
supersymmetric particle (LSP) is not stable and it cannot be a candidate for dark matter. 
The decay of the LSP may produce some interesting signatures at the LHC, which can 
have certain correlations with the neutrino oscillation data. In addition, one can measure 
displaced vertices originating from the LSP decay. 

It should be mentioned here that neutrino masses and bilarge neutrino mixing have also 
been studied in an R-parity violating supersymmetric theory with gauge singlet neutrino 
superfields [O]. However, in that analysis terms of the type v c H\H2 have been dropped 
from the superpotential because of very small coefficient. Analysis has also been carried 
out in the context of the observed baryon asymmetry of the Universe [15|. Finally, one 



should note that the discrete Z% symmetry of the superpotential is spontaneously broken 
in the vacuum. This might lead to cosmological domain wall problem However, the 



solutions to this problem exist [17| and will also work in this case. 
2.2 Soft terms 

Let us now specify the soft-supersymmetry-breaking terms of this model. We will con- 
fine ourselves in the framework of supergravity mediated supersymmetry breaking. The 
Lagrangian £ so ft, containing the soft-supersymmetry-breaking terms is given by 

-£ soft = (mJ)tfQf + (mlfuf % + (m%fdf ~d] + {m\^Lf L* 
+ {mjcfef + m 2 Hl Hf Hf + m 2 H2 Hf H% + (m?-V-W*r,c 



+ Cab 



(A u Y u y j H%QiUj + {A d Y d fHtQ\~d) + (A e Y e f HfL^ + H.c. 
e ab (A u Y v fH b 2 L«vC - e ab (A x \y^H^H b 2 + X -{A K K)^ k v^vl + H.c. 

M3A3A3 + M 2 A 2 A 2 + M1A1A1 + H.c.) . (2.2) 



+ 



1 

~ 2 

In eq,( |2.2| ), the first two lines consist of squared-mass terms of squarks, sleptons (in- 
cluding the gauge singlet sneutrinos vf) and Higgses. The next two lines contain the 
trilinear scalar couplings. Finally, in the last line of eq.(|2.2]), M3,M2, and M\ are Majo- 
rana masses corresponding to SU(3), SU(2) and U(l) gauginos A3, A2, and Ai, respectively. 
The tree-level scalar potential receives the usual D and F term contributions, in addition 
to the terms from £ so ft- 

2.3 The neutral scalar potential and the electroweak symmetry breaking con- 
ditions 

We assume that only the neutral scalar fields develop in general the following vacuum 
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expectation values while minimizing the scalar potential: 



(H° 1 ) = v l , (H° 2 ) = v 2 , (P i ) = u i , 



(2.3) 



The tree level neutral scalar potential looks like 



ijk,.c ci2 
v i v h\ 



(Koutral) = I E Y >i V ] ~ Y + E ' E Y » V ^ ~ X ' V ^ + E K 

i,j i j i i,k 

+ 1 y av^ 2 i 2 + Y \Y Y >^\ 2 + (^j^)E h 2 + n 2 - h 2 ] 2 

i i j i 

+ ™h 2 H 2 +mH>i| 2 + ^(m?)^*^- + ^(m? c )^i/ 



+ 



YiAuYvfwfa - Y(AxX)^!viv 2 + Y \{A KK f k uf V y k + H.c. 



h3 



i,j,k 



(2.4) 



One important thing is to notice that the potential is bounded from below because 
the coefficient of the fourth power of all the eight superfields are positive. We shall further 
assume that all the parameters present in the scalar potential are real. From eq.( |2.4[ ), the 
minimization conditions with respect to Vi,V2, v± are 



2 E^c j + y y ^ 2 + p i v + ma^i + Y^y 1 ^ + (^w = °> ( 2 - 5 ) 



E 



i^C J + Y r c 2 ^ + ^ - 7*&V2 + E(^ y -) iV ^i - (^AA)Vf^ + mf^ = 0, (2.7) 



E A ^ 2 ^ + Ts&vi + - ^Y r i u i + m ^i Vl ~ (AA)Vfw 2 = 0, (2.S 
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where 



(A xX y = E(AWS^ + J2( A « K ) m »W - (Axxyv 1V2 , 

c j = E K * jk ^ + E Y ¥w* - A ^^2, 



,7 = E^W|-(E A x>i' 

/J, = ^2\ l u c i, 

i 

' „2 i „2\ 



7 9 = + ff 2 )> 

i 

^ = EWS-av), 

r^e^E^H 
i 

t# = E^' (2-9) 



In deriving the above equations, it has been assumed that K ij , (A K K)^ k , Y$ , {A v Y v ) l \ 

K) 



(m|o) i - J , (m^yi are all symmetric in i,j,k. 



Note that the Dirac masses for neutrinos are given by = YyV 2 . From present day 
experiments it is well known that neutrino masses are very small. This implies that the 
neutrino Yukawa couplings must also be very small ~ O (10~ 7 ), in order to get correct 
neutrino mass scale using TeV scale seesaw mechanism. This immediately tells us that in 
the limit Y„ J — > 0, eq. ( |2.6D implies that Vi — > 0. So in order to get appropriate neutrino 
mass scale both Y v and Vi have to be small. 

Ignoring the terms of the second order in Y„ J and assuming (vf + vf — v|) ~ (yf — 
= (m~)5 y , we can easily solve eq.(^^) as (using eq. fl2.9Q ) 



f YJ k u k c 3 v 2 - fmYl + (A u Y u yiy 2 \ c j Y^\H x v\ \ 



Note from eq.( 2.10| ), that the left handed sneutrinos can acquire, in general, non- vanishing, 



non-degenerate VEVs even in the limit of zero vacuum expectation values of the gauge 
singlet sneutrinos. However, zero VEVs of all the three gauge singlet sneutrinos is not an 
acceptable solution since in that case no fi-term will be generated. Moreover, one needs 
to ensure that the extremum value of the potential corresponds to the minimum of the 
potential, by studying the second derivatives. 



- 8 - 



3. The scalar sector 



The scalar sector of this model enhances from that of MSSM, because of the choice of 



the superpotential in eq.(2.1) (fourth, fifth and the sixth term). In this case, the neu- 
tral Higgs bosons can mix with both the doublet and gauge-singlet sneutrinos. The CP- 
odd(pseudoscalar) and CP-even(scalar) mass squared matrices are now 8x8, considering 
all three generations of doublet and singlet sneutrinos. Similarly the charged Higgs can mix 
with the charged sleptons and thus the charged scalar mass squared matrix is enhanced to 
8x8. We have considered only the CP-preserving case and hence all the VEVs are chosen 
to be real. The scalar sector of this model has been addressed also in a recent work ||. 
The details of various scalar mass squared matrices are given in appendix [A|. 

For our analytical and numerical calculations in the later part of the paper, we have 
assumed that (m|) lJ = (m|)^, (m? c ) u = (m| c )<P J , and Y^ 3 = 0, Mi / j. We have 
further assumed that K i3k are flavour-diagonal as well as flavour-blind, i.e., K^ k = k if 
i = j = k and zero otherwise. Similarly, we have assumed a flavour-blind coupling A* = A 
for i = 1,2,3. We will see that even with such simplifying assumptions, we can fit the 
global three flavour neutrino data in this model. We will use the following procedure 
for all our subsequent analysis. Using the minimization conditions, we will solve for the 
vacuum expectation values Vi and uf. We will choose the parameters in such a way that 
the values of uf will give an acceptable number for the //-parameter (/i = A]T\i^ c ). As 
a cross check we confirm the existence of two Goldstone bosons in the pseudoscalar and 
charged scalar mass-squared matrices. In addition, we check that all the eigenvalues of the 
scalar, pseudoscalar, and charged scalar mass-squared matrices (apart from the Goldstone 
bosons) should come as positive for a minima. 

Additional constraints on the parameter space can come from the existence of false 
minima. A detailed discussion on this issue has been presented in ref. |J] and the regions 
excluded by the existence of false minima have been shown. One can check from these 
figures that mostly the lower part of the region allowed by the absence of tachyons, are 
excluded by the existence of false minima. In our analysis, we have chosen the parameter 
points in such a way that they should be well above the regions disallowed by the existence 
of false minima. Nevertheless, in the case of gauge-singlet neutrino (v c ) dominated lightest 
neutralino (to be discussed later), the value of k that we have chosen is 0.07 with two 
different values of A, namely, 0.1 and 0.29. In this case, there is a possibility that these 
points might fall into the regions disallowed by the existence of false minima. However, we 
have checked that even if we take the value of k to be higher (0.2 or so), with appropriately 
chosen A, our conclusions do not change much. For such a point in the parameter space, it 
is likely that the existence of false minima can be avoided. A more detailed study on this 
issue is beyond the scope of the present paper. 

Let us also mention here that the sign of the //-term is controlled by the sign of the 
VEV v c (assuming a positive A), which is cotrolled by the signs of AaA and A k k. If ^4aA 
is negative and A k k is positive then the sign of the /i parameter is negative whereas for 
opposite signs of the above quantities, we get a positive sign for the [i parameter. 

The scalar mass-squared matrices (both CP odd and CP even) and the vacuum ex- 
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pectation values vf are not very sensitive to the change in neutrino Yukawa couplings (Y v 
~ O (1CT 7 )) and the corresponding soft parameter A V Y V (~ O (10 -4 ) GeV). On the other 
hand, the values of tan (3 and the coefficients A and k are very important in order to satisfy 
various constraints on the scalar sector mentioned earlier. In fig.|l], we have plotted the 
allowed regions in the (A-/t) plane for tan/3 = 10. 

The values of other parameters are cho- 
sen to be m L = 400 GeV, m&c = 300 GeV, 
y^ 11 = 5.0 x 10~ 7 , Y 22 = 4.0 x 10" 7 , Vj 53 = 
3.0 xlO- 7 , {A v Y v fJ = 1 TeVxy^, (A X X) = 
-1 TeV x A, and (A k k) = 1 TeV x k. The 
upper limit of the value of k is taken to 
be ~ 0.7 because of the constraints com- 
ing from the existence of Landau pole ||. 
With these values of different parameters 
satisfying the constraints in the scalar sec- 
tor, we will go on to calculate the neutrino 
masses and the mixing patterns as well as 
the decays of the lightest neutralino in this 
model as discussed in the next few sections. 

It should be mentioned at this point 
that the radiative corrections to the light Higgs mass, can be significant in some regions of 
the parameter space as discussed in ref.||. It has been shown that the light Higgs mass 
larger than the LEP lower limit of 114 GeV can be obtained with the value of At (trilinear 
coupling in the scalar sector for the stop) within 1-2.4 TeV and when the mixing of the light 
Higgs with the right-handed sneutrino is small. The latter requirement is fulfilled in most of 
the cases that we have considered and in some cases the mixing is slightly larger. However, 
there is always the freedom of choosing the value of At appropriately. Hence, it would be 
fair to say that the experimental limits on the light Higgs boson mass can be satisfied in 
our analysis. The parameter points we have chosen here are sample points with different 
dominant composition of the lightest neutralino. Since we have a large parameter space, 
it is always possible to choose a different parameter point with the same characteristic 
features satisfying all the experimental constraints. 

4. The fermionic sector 
4.1 The neutral fermions 

In this model, because of the breaking of R-parity, two neutral gauginos, B°(= and 
Wg (= — iA 3 ,) and two neutral higgsinos (H® and H®) are now mixed with the neutrinos 
(both Vi and vf). As can be seen from the superpotential (fourth and fifth term in eq. 
d2.1|) ), the fermionic partners of vf and z/j mix with the neutral higgsinos ||, The neutral 
gauginos are mixed with the left-handed neutrinos through the vacuum expectation values 
of the doublet sneutrinos. Mass matrices for the neutral and charged fermion sectors, 




U01 0.05 0.10 0.15 0.20 0.25 0.30 



Figure 1: Allowed regions in (A-k) plane which 
satisfy various constraints on the scalar sector, 
for tan/3 = 10. A and k were allowed to vary from 
0.005 to 0.50 and 0.005 to 0.70, respectively. 
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involving all three generations of neutrinos (both doublet and singlet) and charged leptons, 
have been addressed also in ||. 

In the weak interaction basis defined by 



B°, Wl Hi Hi u c e , vl u e , u„, u T ) , 



the neutral fermion mass term in the Lagrangian is of the form 



^neutral 



^ oT M n ^ + R.c., 



(4.1) 



(4.2) 



where A4 n includes all three generations of doublet and gauge singlet neutrinos and thus 
it is a 10 x 10 matrix. The massless neutrinos become massive due to this mixing with the 
neutralinos and the gauge singlet neutrinos. The three lightest eigenvalues of this 10 x 10 
neutralino mass matrix correspond to the three light physical neutrinos and their masses 
have to be very small in order to satisfy the experimental data on massive neutrinos. The 
matrix A4 n can be written in the following manner 



M r 



M 7x7 "4x7 
"1-3x7 3X 3 



(4.3) 



where 



and 



/ M 1 




M 



7x7 





M 2 



-f 2 v 2 



\ 




V^ 1 T2 Vl 
91 



92 



-\ e V2 —X^V2 —X T V2 



V2 






l>2 



92 






v 2 



-\ e V 2 P e 



2uT 2uT 2u e J 



-\»V2 2u 



fie 



2u\ 



2nr 



V o 



2ul e 2uS 2uV 



m 3x7 



X T v 2 p 

f~fae %Ve T% Y u ee V 2 Y^V 2 Y?V2\ 

fa, fan r£ Yrv 2 Yrv 2 Y^ T v 2 



(4.4) 



(4.5) 



Note that the top- left 4x4 block of the matrix M? X 7 is the usual neutralino mass 
matrix of the MSSM. The bottom right 3x3 block is the Majorana mass matrix of gauge 
singlet neutrinos, which will be taken as diagonal in our subsequent analysis. The entries 
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of Mjx7 are i n general of the order of the electroweak scale and the entries of m^ x j are 
much smaller (~ 0(1O -5 GeV)). Hence, the matrix ( |4,3[ ) has a seesaw structure that will 
give rise to three very light eigenvalues corresponding to three light neutrinos. The correct 
neutrino mass scale of ~ 10~ 2 eV can easily be obtained with such a structure of the 
10 x 10 neutralino mass matrix. In this work our focus would be to see if one can obtain 
the correct mass-squared differences and the mixing pattern for the light neutrinos even if 
we consider flavour diagonal neutrino Yukawa couplings in eq.( [4.5| ) (i.e. with a diagonal 
Dirac neutrino mass matrix). This makes the analysis simpler with a reduced number of 
parameters and makes the model more predictive. As we will show later, it is possible to 
find out the correct mixing pattern and the mass hierarchies (both normal and inverted) 
among the light neutrinos in such a situation, at the tree level. 

In order to obtain the physical neutralino states, one needs to diagonalize the 10 x 10 
matrix M n . As in the case of MSSM, the symmetric mass matrix A4 n can be diagonalized 
with one unitary matrix N. The mass eigenstates are defined by 

X° = iVJ>;, »,j = l,...,10, (4.6) 

where the 10 x 10 unitary matrix TV" satisfies 

N*M n N- 1 = M° D , (4.7) 

with the diagonal neutralino mass matrix denoted as The matrix N may be chosen in 

such a way that elements of are real and non-negative. In our analysis we will assume 
that all the entries in the 10 x 10 neutralino mass matrix M n are real. Seven eigenstates 
of this matrix are heavy, i.e. of the order of the electroweak scale. Out of these seven 
states, there are four states which are usually very similar to the MSSM neutralinos. The 
remaining three states are mostly dominated by vf. It is, in general, very difficult to predict 
the nature of the lightest of these seven states since that depends on several unknown 
parameters. In our analysis of the decays of the lightest supersymmetric particle (LSP), 
we will concentrate on three different possibilities: (i) lightest state is dominated by the 
bino component, (ii) higgsino dominated lightest state, and (iii) gauge singlet neutrinos vf 
form the lightest state. The last possibility is very interesting since in this case we have the 
opportunity to produce the gauge singlet neutrinos at the LHC and study their properties 
through the R-parity violating decay modes. This way one has a direct probe to the seesaw 
scale at the LHC. 

4.2 The charged fermions 

In the charged fermion sector, the charged gauginos and charged higgsinos mix with the 
charged leptons because of the presence of the effective bilinear RPV parameters e l = 
Y^j Y^Vj and the sneutrino VEVs z/j. This is similar to the case of MSSM with bilinear 
RPV with the parameter fi defined as \i = ^ AV?. Since we want to calculate the decays of 
the lightest neutralino, we also need to know the mass eigenvalues and the mixing matrices 
in the charged fermion sector. Because of this reason, we discuss the chargino mass matrix 
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in some details. In the weak interaction basis defined by 



VT 



-T 



( e R , t r ), 



The charged fermion mass term in the Lagrangian is of the form 



pmass 
^charged 



T 1 T T T 



( 5 



T 

x5 m 5x5 



y m5 X 5 05x5 



(4.9) 



Here we have included all three generations of charged leptons and assumed that the 
charged lepton Yukawa couplings are in the diagonal form. Also, — iX^ are the two- 
component charged Wino fields and and are the two-component charged higgsino 
fields. The matrix 7775x5 is given by 



/ M 2 92^2 











-Ytve -Y^y -Y^v T 



"75x5 



g 2 v e ~r% Y* e vi 



(4.10) 



92V u -r£ Yt^vi 



\92Vr ~r T c XTn J 

The charged fermion masses are obtained by applying a bi-unitary transformation such 
that 

U*m^V- 1 = M% (4.11) 

where U* and V are two unitary matrices and Ai D is the diagonal matrix with non-negative 
entries corresponding to the physical fermion masses. The two-component mass eigenstates 
are defined by 

xr = Uij*J, i,j = l,...,5. (4.12) 

Nevertheless, we notice that the 13, 14, and 15 elements of the chargino mass matrix (eq. 
( 4.10D ) are vanishing and given the orders of magnitude of various parameters, we also see 
that the values of the other off-diagonal entries (except for 12 and 21 elements) are very 
small. This indicates that the physical charged lepton eigenstates will have a very small 
admixture of charged higgsino and charged gaugino states. So we can very well assume 
(also verified numerically) that this mixing has very little effect on the mass eigenstates 
of the charged leptons. Thus, while writing down the neutrino mixing matrix, it will be 
justified to assume that one is working in the basis where the charged lepton mass matrix 
is already in the diagonal form. 
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5. Neutrinos 



5.1 Seesaw masses 

When lepton-number violation is allowed, the effective light neutrino mass matrix arising 
via the seesaw mechanism is in general given by 

M u = -m 3x7 M^ 7 mJ x7 . (5.1) 

As discussed in the previous section, 7713x7 is the so-called Dirac neutrino mass matrix 
and Mjxj is the matrix for the heavy states and contains AL = 2 mass terms for right chiral 
neutrinos. In order to find out the neutrino mass-squared differences and mixing angles, 
one must diagonalize the matrix M u to find out the eigenvalues and the eigenvectors. 
Before discussing the detailed numerical results, let us try to understand the characteristic 
features of this neutrino mass matrix analytically. Let us note that the neutrino mass- 
squared differences indicate three possible scenarios for the light neutrino mass spectrum. 
They are (i) Normal hierarchy corresponding to mi ~ 771-2 ~ \J ^ m 2i 1 m 3 ~ \/i^^3iT> ( n ) 
Inverted hierarchy : mi ~ 7772 ~ ylAmgJ , 7773 <C yfAmf]] and (iii) Degenerate masses : 
mi ~ m 2 ~ m 3 3> -\/|^^r7^3xl , where mi, 7772, and 7773 are the three light neutrino mass 
eigenvalues. 

5.1.1 Analytical results 

Let us make a few simplifications in order to diagonalize the effective light neutrino mass 
matrix. One should note that the neutrino mass matrix involves the vacuum expecta- 
tion values for the doublet and the gauge singlet sneutrinos. Hence we also make some 
simplifying assumptions for the parameters appearing in the scalar sector. Some of these 
assumptions have already been mentioned in the scalar sector but here we repeat them for 
the convenience of the reader. We have defined 

K tjk = Kj ^ i = j = an( j zero otherwise, 
{A K K)^ k = (A k k), if i = j = k, and zero otherwise, 
Y?=0, iii^j, 
(A v Y v f = 0, if i + j, 
A 1 = A 2 = A 3 = A, 

(A x xy = (A X X) 2 = (A A A) 3 = (A\X), 
{m\) ij = (m 2 L )5 ij , 

(mlf = (ml)^, (5.2) 

where i, j, k = e, /j, r in the flavour basis. 

With these assumptions one can solve for from the minimization equations for the 
gauge singlet sneutrinos (eq. (|2.5| )) and the result is v\ = = ^3 = v c . This can be 
understood if we neglect the terms proportional to Y u 2 , Y v v, {A y Y v )v in the minimization 
equations and assume that X,k are ~ 0(1) couplings. 
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Now let us look at the effective left chiral neutrino mass matrix (eq,( |5.1|) ) in a little 
more details. Because of the smallness of Y v and Ui, one can neglect terms containing 
Y 2 v 2 and Y v z v. This way one obtains an approximate analytical expression for the 3x3 
neutrino mass matrix. 



where 



2Av c 



3 A 



/ b 2 e b e bn b e b T \ 



b e b^ b 2 ^ b^b 7 



+ 



6kv c 



\ b e b T b^ b 2 T J 



( -la 2 a e aa a e a T \ 



a e a^ -2a M a M a T 



\ a e a T a^a T -2a 2 j 



( 



2A M v- vii ( 2\AB\ 



c e c e c^ c e c T ^ 



\ c e c T C^C T C T J 



A = \ z {v{ + v'if + A\kv c2 Vi v 2 - AMXAfji, /x = 3Ai/ c , 
A = (kv c2 + \vxvi), 
1 

M ~ 



Mi M 2 



B = vi{vl + vl) -2MHV2, 

<H = Y?V2, h = (Y*vi + 3\ui), a = ^ 



(5.3) 



(5.4) 



with k = e, fi, r. 

One can rewrite eq. (|5 



Ml 



in a compact form as follows 



1 



2Av c 

— —bib,-, H aiaJl — 35{ 

3A 3 §kv c 3K 



2 Vy y . 

3 A ^ l 



2XAB\ 



(5.5) 



Let us note that the smallness of the left chiral sneutrino VEVs (z/j << t>i,t>2) allows 
us to use ra| « \{g\ + fli)( u i + u l) an d tan /? ~ 

The coefficients of the first term in eq. (|5.3| ) (or in eq.(|5,5|)) is of the order of A whereas 
the coefficient of the second term is ^ iq^jj where m is the electroweak (or supersymmetry 
breaking) scale and we have assumed that the relevant mass scales are at fh and k is an 
order one coupling. The value of the coupling A (determines the value of the \i parameter), 
which satisfies the neutrino data as well as the constraints in the scalar sector, is taken to 
be of the order of 10 _1 . On the other hand, the coefficient of the third term is ~ -~-(^). 
Since j£ ~ 10 _6 -10 , there is an extra suppression factor in the elements of third term 



in eq.( |5.3|) , compared to the first two terms. In addition, b 2 ~ a 2 ~ c? with a slightly 
larger value of af and cf compared to bf in most cases. Hence, one can neglect the third 



term of eq.(5.3) in comparison to the first two terms. However, in our numerical analysis 



- 15 - 



(discussed later) we have kept all the terms in eq. (|5.3| ) and checked that the presence of 
the third term changes the result in an insignificant manner. 

Before going on to find out the expressions for the eigenvalues and the eigenvectors of 
the effective neutrino mass matrix (eq.( |5.3| )), let us highlight a few limiting cases which 
give us some insight regarding the behaviour of the neutrino mass matrix. Neglecting the 
third term one can rewrite eq,Q5.5|) in the following manner 



J bnu c 



Y?Y'i(l-mj) 



1 

2M 



UiUj + 



vi v c (Y™ Vj + YP Vi ) Y* Yi 3 v\v 



2,,c2 



+ 



— — — kv c1 sin 23 H 

2MAn V 2 



(5.6) 



Here we have used V2 = v sin/3, v\ = v cos 3, and [i = 3Xu c . 
In the limit v c — > oo and v — > 0, eq. (EM) reduces to 



' 2M' 



(5.7) 



which is the first part of the second term of eq.(|5.6|). In this case the elements of the 
neutrino mass matrix are bilinears in the left-handed sneutrino VEVs and they appear 
due to a seesaw effect involving the gauginos. This is called the "gaugino seesaw" effect 
and neutrino mass generation through this effect is a characteristic feature of the bilinear 
R— parity violating model. This effect is present in this model because we have seen earlier 
that the effective bilinear R-parity violating terms are generated in the scalar potential as 
well as in the superpotential through the vacuum expectation values of the gauge singlet 
sneutrinos. Note that the gaugino seesaw effect can generate mass for only one doublet 
neutrino. 

In the limit M — > oo, eq. (15.61) reduces to 



M v - 



Qkv' 



: Y?YV(l-36ij), 



(5.8) 



which corresponds to the ordinary seesaw effect between the left handed and gauge singlet 
neutrinos. Remember that the effective Majorana masses for the gauge singlet neutrinos are 
given by Mr = 2kv c . The ordinary seesaw effect can generate, in general, masses for all the 
three neutrinos. Thus depending on the magnitudes and the hierarchies of various diagonal 
neutrino Yukawa couplings Y™, one can generate normal or inverted hierarchy of neutrino 
masses (combining with the "gaugino seesaw" effect) corresponding to atmospheric and 
solar mass squared differences, as discussed earlier. In this model it is difficult to obtain 
a degenerate neutrino spectrum and we do not consider this possibility in our subsequent 
analysis. 

Now let us try to find out the approximate analytical expressions for the eigenvalues 
and eigenvectors of the effective light neutrino mass matrix using perturbation theory. 
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Neglecting the third term in eq,(|5.3|), the neutrino mass matrix looks like 



M u = B 



( bl 



b e b^ b e b T \ 



bebfj, 



bab T 



+ A 



\ b e b T b^ b T j 



I -2a 2 o e a M a e a T \ 



a e a^ -2a M a^a T 



\ a e a T a M a r -2a 2 / 



(5.9) 



where A = and B = §4|~- As we have argued above, the first matrix in eq.(|5.9|) can 
be considered as the unperturbed one and the second matrix can be treated as a perturba- 
tion over the first one because of the presence of the smaller coefficient A. The eigenvalues 
of the unperturbed matrix are (0,0, /3(6 2 + b 2 + b 2 )) and the corresponding eigenvectors 

1 I . With the order of magnitudes of various 



arc 



-fc 01 



£ 1 



r 



b T b T 



parameters discussed above, the only non-zero eigenvalue determines the atmospheric neu- 
trino mass scale corresponding to the normal hierarchical mass pattern for neutrinos. In 
order to generate the solar neutrino mass scale one must turn on the perturbation. In this 
case one should use the unperturbed eigenvectors to get the corrections to the eigenvalues 
due to perturbation. However, since two of the eigenvalues are zero, one needs to apply 
the degenerate perturbation theory to evaluate the correction to the eigenvalues. To do 
this, first we construct a complete set of orthonormal eigenvectors using Gram-Schmidt 
orthogonalization procedure. The set of orthonormal eigenvectors obtained in this case are 



Y2 



vw+w 

/ be 



br 

V 1 



be 


V 1 / 



1 

\ ~bj+bj 



(5.10) 



where 



bl + bl + bl. 



(5.11) 



Using degenerate perturbation theory for this set of orthonormal eigenvectors, the 
modified eigenvalues m'± and 777,3 are obtained as 



m'± 



777o 



A_ 



sn 2 h 



n afe ±^[-30 2 (£ afe ) 2 + (n a 

^|(E a ^) 2 -3A afc |, 



(5.12) 
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where 



A a fe = / aiajbjbj, 

i<j 

n a f, = ^(aibj + ajbi) 2 - A ab , 

i<j 

Safe = ^2 aiajb k . (5.13) 

As one can see from eq.( [5.12| ), the corrections to the eigenvalues are proportional to the 
coefficient A appearing in eq. ( [5.9D . This is the effect of the ordinary seesaw. Let us note in 
passing that this effect is absent if only one generation of left chiral neutrino is considered, 
whereas for two and three generations of left chiral neutrino the ordinary seesaw effect 
exists. This can be understood from the most general calculation involving n-generations 
of left chiral neutrinos, where the coefficients of A pick up an extra factor (n — 1). 



5.1.2 Numerical Results 

In order to get some idea about the numbers of the mass eigenvalues and to make com- 
parisons between the full numerical results and the results using approximate analytical 
expression, we look at a sample point in the parameter space. As mentioned earlier, 
we are considering only the normal hierarchical pattern of neutrino masses. The set of 
parameters are M x = 325 GeV,M 2 = 650 GeV, A = 0.06, k = 0.65, A\\(— A k k) = 
-ITeV x A(k) and tan/? = 10. 

The choices of diagonal neutrino Yukawa couplings (Y™) and corresponding soft param- 
eters (A u Y u y i are very crucial and we take, for this particular calculation, Y£ e = 4.57x 10 -7 , 

= 6.37 xl0~ 7 , YJ T = 1.80 xl0~ 7 , (A u Y u ) ee = 1.57xl0~ 4 GeV, (A U Y U )^ = 4.70 xl0~ 4 
GeV, {A U Y V ) TT = 3.95 x 10~ 4 GeV. Soft masses of left handed and right handed sleptons 
are chosen to be 400 GeV and 300 GeV, respectively. Later on we will show the allowed 
regions in the Y v planes which satisfy the experimental data on neutrino masses and mix- 
ing. For these choices of various parameters, the derived left-handed sneutrino VEVs are 
v e ~ 10~ 5 GeV, = 1.515 x 10~ 4 GeV, v T = 2.133 x 10" 4 GeV and right-handed sneutrino 
VEVs are v c = —588.74 GeV. With this set of values the masses of three neutrinos have 



been found out by direct diagonalization of the matrix obtained using (5.1) and also from 



the approximate analytical expression using ( 5.12 ). It has been observed that even with 



several simple assumptions (eq.( |5.2| )), all three generations of left chiral neutrinos acquire 
non- vanishing, non-degenerate masses at the tree-level. The comparison of the results as 
obtained from ( |5.1| ) and from ( |5.12| ) are given in Table |l[ One can see that these values 
are within the 3a limits shown in eqs.([l.l|)-( |1.2| ). However, it should be mentioned that 
if A, k are much less than ~ Oil) and Y u s are much larger than the ones we have consid- 
ered above, the approximate analytical expression does not produce the correct results for 
the eigenvalues and the eigenvectors. Obviously, when the neutrino Yukawa couplings are 
larger one cannot consider the second term in eq. (|5.3| ) as a perturbation to the first term. 
In our numerical analysis for obtaining the allowed region of parameter space which satisfy 
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(eV) (xlO 3 ) 


Am^^eV 2 ) 
(xlO 5 ) 


Amf^eV 2 ) 
(xlO 3 ) 


mi 


m 2 


m 3 


eq.© 


9.970 


4.169 


48.23 


8.203 


2.307 


eq.(|5.12D 


9.468 


4.168 


47.71 


7.228 


2.187 



Table 1: Absolute values of the neutrino masses and the mass-squared differences for a sample point 
of the parameter space discussed in the text. Results for full numerical analysis have been obtained 
using eq.(pTTj). Approximate analytical expressions of eq.(5.12) have been used for comparison. 



the neutrino data, we have done a full numerical analysis without using the approximate 
formula. 

The numerical values of the solar and atmospheric mass squared differences Am^ 
and Anijj have also been shown in Table |l| and the results show good agreement. The 
numerical calculations have been performed with the help of a code developed by us using 



Mathematica [18]. In our numerical calculations, we have taken for (i) normal hierarchy: 
iri2\max < 1-0 x 10 ^ GeV and (ii) inverted hierarchy: Tn^\ max < 1.0 x 10- 11 GeV. 



5.2 Neutrino mixing 

The left chiral light neutrinos form a 3 x 3 mass matrix in the flavour basis. The unitary 
matrix which diagonalizes this mass matrix can be parameterized as follows |19|], provided 
that the charged lepton mass matrix is already in the diagonal form 



C12C13 



-S12C23 - Cl2S23Sl3e" 



\ S12S23 - Cl2C23«13e 



•S12C13 si3e" 



C12C23 - Sl2S23Sl3e S23C13 



-Cl2«23 - Sl2C23Sl3e * C23C13 / 



(5.14) 



where c 



COS I 



sin(9jj, and i,j runs from 1 to 3. Various neutrino oscillation 



experiments indicate that #12 ~ 34°, #23 ~ 45°, and #13 < 13° [20, 21]. This pattern is 



known as bilarge mixing. In order to understand the consequences of such mixing in the 
zeroth order, one can approximately take #23 = 45°, sin #12 = ^= and #13 as 0°, something 



known as tribimaximal structure [22]. Then the unitary matrix turns out to be 



U 



3x3 



1 



1 1 

%/3 V2 



(5.15) 



V — -— — / 



Given the three mass eigenvalues mi,m2,m3, it is possible to use the matrix U. v to obtain 
the mass matrix in the flavour basis as follows, 
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W M u U h 



M. 



diag ' 



(5.16) 



where 



^diag 



(mi \ 



m 2 



V m 3 / 



(5.17) 



We will numerically diagonalize the neutrino mass matrix M v obtained in eq.([|[l]) and also 
use the approximate analytical method to find out the neutrino mixing matrix U v . We 
will also compare the results obtained using these two methods. However, when we will 
scan the parameter space to find out the allowed regions where the neutrino experimental 
data are satisfied, we shall use the full numerical procedure. The advantage of having 
the approximate analytical expression is that it can give us some insight regarding the 
conditions on the model parameters for which the bilarge mixing is obtained. We can 
verify these predictions numerically in some regions of the parameter space. We will try 
to find out the regions in the models parameters where the numbers in eqs.([Ql)- (|l.3|) are 
reproduced. 



5.2.1 Analytical results 

One can construct the neutrino mixing matrix analytically using the degenerate perturba- 
tion theory. With the set of orthonormal eigenvectors in eq. ( |5.10| ) and the eigenvalues in 
eq.( 5.12j ), it is possible to write down the eigenvectors of ( |5.9| ) in the following form 



(y 



l)3xl = «12/1 + «22/2, 



(5.18) 



(^2)3x1 

(^3)3x1 



a'iyi + a' 2 y2, 



2/3, 



(5.19) 
(5.20) 



where a±, 02, a[, a' 2 , are calculated using degenerate perturbation theory and their ana- 
lytical expressions are given by 



Oil 



h 



12 



(5.21) 



a 2 




(5.22) 



(5.23) 
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/ 

Q 2 



/ill — m '- 



\Jhl 2 + (h n -m'_f 
Here rn'i, m'_ are given by eq.(5.12) and tin, h\2 ar e given by 

2A {a 2 b 2 + a e a T b e b T + a 2 ^ 2 -) 



(5.24) 



and 



where 



H2 



A [a^(a T b e - a e b T )b\ - b^ [2b e b T a 2 _ + a e a T b 2 _)] 



n h bl 



(bl±b 2 T ), 



(5.25) 
(5.26) 



(a 



(5.27) 



and Qf, has been defined in eq.( ^.ll ). 



The neutrino mixing matrix U u can be constructed using these eigenvectors in eqs.( |5.1i 



( 5.20 ) and it looks like 



yi y2 y 3 



3x3 



(5.28) 



Looking at the expressions for the eigenvectors, one can immediately draw a few conclusions 
regarding the behaviours of the neutrino mixing angles with the model parameters. For 
example, the (13) mixing angle #13 is given by 

hi 



sin 0i3 



b 2 + b 2 +b 2 T - 



If we want the (13) mixing angle to be small then one must have b 2 <C 
other hand, the (23) mixing angle #23 is given by 



sin t/ 2 3 



bi+br 



(5.29) 
On the 

(5.30) 



So, if the (23) mixing is maximal then one would expect b 2 l = b 2 . The solar mixing angle 
#12 is approximately given by 



be 

sin 2 #i2 ~ 1 — (at'i + a 2 7~) 

T 



(5.31) 



where and a' 2 are given by eqs,( 5.23| ) and ( 5.24] ), respectively. In order to have #12 ~ 
35°, the square root of the second term on the right hand side of eq.( |5.31 ) should be 
approximately 0.8. In the next sub-section we discuss the patterns of neutrino mixing in 
this model numerically, and show the allowed regions of the parameter space where the 
neutrino experimental data are satisfied. 
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5.2.2 Numerical Results 

Let us first calculate the neutrino mixing 
angles for the parameter point discussed in 
Table |. As we have discussed earlier, this 
parameter point generates the normal hi- 
erarchical pattern of neutrino masses. In 
Table ^, the three mixing angles are shown 
and they have been evaluated using the di- 
rect numerical calculation in eq.(5.1) as well 
as using the approximate analytical expres- 
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of other parameters are described in the text. The 
lightest neutralino (LN) is either bino or higgsino 
dominated. 



sions in eq.( 5.28 ). We want to emphasize 
once again that the approximate formulae 
have been used just to get some idea about 
the behaviours of the neutrino masses and 
mixing angles with various model parame- 
ters. This way we can identify the relevant 

parameters which crucially control the neutrino masses and mixing angles in different re- 
gions of the parameter space. However, these formulae are not valid everywhere in the 
allowed parameter space and in all the plots shown in this paper we have used full numer- 
ical calculation using eq.(|5.1|). 



mixing angles in degree 


Using (|5T|) 


Using (|5.280 


#12 


36.438 


37.287 


#13 


9.424 


6.428 


#23 


38.217 


42.675 



Table 2: Neutrino mixing angles for a sample parameter point discussed in Table |]. Results are 
shown using eq.(5.1) and eq.(5.2E). 



We have taken suitable values for A and k in such a way that they fall in the region 
allowed by the constraints in the scalar sector (similar to the region shown in fig. |l|). We can 
see that for this choice of the parameter space, numerical and approximate analytical results 
give quite good agreement. Naturally, one would be interested to check the predictions 
made in eqs. ( PT29D , flOcD , and flOU ) over a wide region in the parameter space and see 
the deviations from the full numerical calculations. This has been shown in fig. [|, where 
we have plotted the value of sin 2 023 as a function of the ratio 6 2 /& 2 . 

We can see from this figure that for 6 2 = 6 2 , the value of sin 2 #23 varies in the range 0.41 
- 0.44, which corresponds to #23 between 40° and 42°. On the other hand, eq,( |5.30 ) tells 
that for 6 2 = 6 2 , sin 2 #23 = 0.5. So we see that in this case the result from the numerical 
calculation is reasonably close to the prediction from the approximate analytical formula. 
The choices for various parameters are given below. For the gaugino mass parameters 
Mi and M2, we take two different sets of values which give us either a bino dominated 
lightest neutralino or a higgsino dominated lightest neutralino. In order to have a bino 
dominated lightest neutralino, our choices are Mi = 110 GeV,M2 = 220 GeV, and for 
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a higgsino dominated lightest neutralino we take M\ = 325 GeV, M2 = 650 GeV. The 
same choices will be made for the gaugino mass parameters when we discuss the decays of 
the lightest neutralino in Sec. VI. Our choice of the ratio of the gaugino mass parameters 
at the electroweak scale is motivated by the assumption of universal gaugino mass at the 
grand unified theory scale. The value of k is taken as 0.65 which satisfies the constraints 
from the scalar sector. We have taken two different values of A corresponding to a bino or 
a higgsino dominated lightest neutralino. For the bino dominated case A = 0.13, and for 
the higgsino dominated case A = 0.06. The corresponding values for A\\ = —A x ITeV 
and A k k = k x ITeV. The three diagonal neutrino Yukawa couplings (Y™) vary randomly 
in different ranges 

3.55 x 10~ 7 < Vj 1 < 5.45 x 10~ 7 
5.55 x 10~ 7 < Y 22 < 6.65 x 10~ 7 

1.45 x 10~ 7 < Y™ < 3.35 x 10~ 7 . (5.32) 

The corresponding soft parameters (A™) also vary randomly in different ranges such that 
the parameters (A u Y u ) n effectively vary as follows 

1.25 x 10~ 4 < {A U Y U ) 11 < 1.95 x 10~ 4 
3.45 x 10" 4 < {A V Y V ) 22 < 4.95 x 10~ 4 

2.35 x 10~ 4 < (A u Y u f 3 < 4.20 x 10" 4 . (5.33) 

The allowed regions in the A — k plane are not very sensitive to the values of Y v and 
A V Y V due to their smallness. Hence we choose them different for different cases, in order 
to accommodate the three flavour global neutrino data. 
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a 2Tr) 2 - One can see that as (a^ +a2^) 2 ~~ * 
0.50, sin 2 9i2 tends to 0.50, as predicted by 
the analytical formula. Values of other pa- 
rameters are the same as in fig. 2. 



Note that in some parts of these ranges we have considered a bino dominated lightest 
neutralino and in some parts we have taken a higgsino dominated lightest neutralino with 
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some overlapping regions. The values of other parameters are chosen to be = 400 GeV, 
m^c = 300 GeV and tan/3 = 10. We have assumed that the phase 5 appearing in the 
mixing matrix ( |5.14 ) is zero. One important thing to notice is that even with a flavour 
diagonal structure of the neutrino Yukawa couplings Y v , one can obtain the required two 
large mixing angles for the neutrinos. The variations of other two mixing angles with the 
relevant parameters are shown in Figs. ^ and 

In fig|| we have shown the regions in the various Y v planes satisfying the three flavour 
global neutrino data. The values of other parameters are as in fig.|2| for the case where 
the lightest neutralino is bino dominated. We can see from these figures that the allowed 
values of Y u s show a mild hierarchy such that Y 22 > Y^ 11 > Y^ 3 . 

Similar studies have been performed for the inverted hierarchical case and the allowed 
region shows that the magnitudes of the neutrino Yukawa couplings are larger compared 
to the case of normal hierarchical scheme of the neutrino masses with a different hierarchy 
among the Y u s themselves (Y^ 11 > Y 22 > Yj? 3 ). In this case sin 2 6±2 shows an increasing 
behaviour with the ratio 6 2 /& 2 , similar to the one shown by sin 2 623 with b^/b 2 in the normal 
hierarchical scenario. On the other hand, sin 2 623 shows a decreasing behaviour with b^/b 2 . 
In all these cases, the solar and atmospheric mass-squared differences are within the 3a 
limits. 



The case of v c dominated lightest neutralino has also been studied and it shows a very 
interesting and different behaviour compared to the bino and higgsino dominated cases. 
In this case, the dominant contribution in the neutrino mass matrix (eq.( |5.5D ) comes from 
the term proportional to a%aj. The terms proportional to bibj should be considered as a 
perturbation. Hence, in the normal hierarchical scenario of neutrino masses, one would 
expect that sin 2 623 is proportional to a 2 /a 2 . This is exactly what we see in fig.||. 

Note that for a 2 = a 2 , the mixing becomes maximal. On the other hand, the solar 
mixing angle is controlled mostly by the quantity fe 2 /^ and shows an increasing behaviour 
with this ratio. In the case of inverted hierarchical scenario of neutrino masses, sin 2 623 
shows a decreasing behaviour with the ratio b^/b 2 whereas sin 2 612 shows an increasing 
pattern with 6 2 /6 2 . However, we do not show these plots here. 



6. Decays of the lightest neutralino 

Let us now look at some decay processes which can be considered as the typical consequence 
of this model. It is obvious that because of the R-parity violation there will be no stable 
lightest supersymmetric particle (LSP) present in this model. Here we consider the case 
where the lightest neutralino (x? in our notation to be described below) is the LSP (or 
NLSP in some cases) and study its decay pattern in the R-parity violating channels. In 
particular, we will consider the case where m^o > m^±, so that the three-body decays 
are less important compared to the two-body decays x.7 ~ * % + u e,n,ri X7 ~> + eT i 
X® — ► + fi^, and X7 ~^ W + T=F - The required Feynman rules for the computation 
of the decay of the lightest neutralino are given in the Appendix [B]. Let us also remark 
that the lightest neutralino LSP can also decay to h + u, if it is kinematically allowed, 
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Figure 5: Plots for normal hierarchical scheme of neutrino mass in Y„ ~ Y 22 , Y^ 1 — Y 33 and 
Y 22 — Y 33 plane when the lightest neutralino (LN) is bino dominated. 

where h is the MSSM-like lightest Higgs boson. However, for our illustration purposes 
we have considered the mass of the lightest neutralino in such a way that this decay is 
either kinematically forbidden or very much suppressed (assuming a lower bound on the 
mass of h to be 114 GeV). Even if this decay branching ratio is slightly larger, it is usually 
smaller than the branching ratios in the (£{ + W) channel. Hence, this will not affect our 
conclusions regarding the ratios of branching ratios in the charged lepton channel (£i + W), 
to be discussed later. The lightest neutralino decay x? ~~ * v + where v c is the scalar 
partner of the gauge singlet neutrino v c , is always very suppressed. We will discuss more 
on this when we consider a u c dominated lightest neutralino. 
Consider the following decay process 

Xi^Xj + V, (6.1) 

where Xi(j) ^ s either a neutralino or chargino, with mass and V is the gauge boson 

which is either or Z, with mass m v . The masses m; L and rrij are positive. 
The decay width for this process in eq. (|6.1D is given by [23, 24] 



T (xi — Xi + V)= J 2}C T 2 x{{G 2 L + G 2 R )T-GlG R g}, 
32 irmfm^y 
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where J-, Q are functions of , nrij , m v and given by 

T{mi, m,j,m v ) = K + 3 m 2 , (mf + m 2 — m 2 ,) , 
eiejmimjml, 

(6.2) 



Q(rn,i,mj,m v ) = 12 eiejmimjml 



with £j(j) carrying the actual signs (±1) of the neutralino masses. The chargino masses 
must be positive. The kinematical factor KL is given by 

JC(m 2 , m 2 , m 2 ,) = (mf + m 2 — m 2 ,) 2 — 4 mfm 2 . (6-3) 



In order to derive eq. ( |6.2[ ) , we have used 
the relation m\y = m 2 z cos 2 #w and since 
(i/i) « v\,V2, some of the MSSM relations 
still hold good. The factors Gl, Gr are given 
here for some possible decay modes 

For Xi^XjZ 
Gl = Oj^, Gr = OjP, 
For Xi^XjW- 
Gl = Ofj, Gr = Ofj, 

(6.4) 

where ofj and 0"f are given by eq.( |B.5| ) 
and eq.Q). 

Now consider the following decays 
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Xln — > Z + u, 
Xln^W ± + £T, 

(6.5) 

where x°ln stands for lightest neutralino and £ = e, /i, r. At this stage let us discuss our 
notation and convention for these decays. The neutralino mass matrix is a 10x10 mass 
matrix which includes the left handed as well as the gauge-singlet neutrinos. If the mass 
eigenvalues of this matrix are arranged in the descending order then the three lightest 
eigenvalues of this 10x10 neutralino mass matrix would correspond to the three light 
neutrinos. Out of the remaining seven heavy eigenvalues, the lightest one is denoted as 
the lightest neutralino. Thus, in our notation X7 is the lightest neutralino and x^ +7 , where 
j = 1, 2, 3 correspond to the three light neutrinos. Similarly, for the chargino masses, x^j_ 2 
(I = 1, 2, 3) corresponds to the charged leptons e, fx, r. 
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(6.6) 



(6.7) 



So for x\n Z + u, which is also x 7 - *• Z + Xj+7 U = 1; 2, 3), one gets from eq. (|6~< 
and eq.( |B.5| ) 

Gl = ~N j+7)3 N^ + \n 3+7A N*, - ^N j+7tk+r Nf >k+7l 
where j, k = 1, 2, 3 and this in turn modifies eq.( |6.2j ) as 

r (x°r - z + x° j+ r) = „ 9 g2]C T s x {2 + g£ a} , 

J 32 irmlnin?,, L J 

with mj = m^o, = ~ and m v = mz- 

Let us now consider the other decay which is x°ln ~~ * ^ ± + ^ or equivalently x 7 
W ± + xf(j = 3,4,5). 

For the process x° — > + x+ 

r (x? - + if) = JXni, x < < Gi + °«> ' - Gla " e > • 

x" w 

G L = N 72 V^ - -^N U V* 2 , 

G R = Nyj jX + -^7*3^2 + -±=Nl k+7 U jjk+2 , 

(k = 1, 2, 3), (6.8) 

where eq.(6.4) and eq.( [B.8D has been used. The process x 7 — * W + + Xj is obtained by 
charge conjugation of the process in eq.(|6.8j). 

6.1 Correlation between the lightest neutralino decays and neutrino mixing 
angles 

Correlations between the lightest neutralino decays and neutrino mixing angles will depend 
on the nature of the lightest neutralino as well as on the mass hierarchies of the neutrinos, 
i.e. whether we have a normal hierarchical pattern of neutrino masses or an inverted 
one. In this section we look into these possibilities in details and consider three different 
cases for the dominant component of the lightest neutralino. We consider that the lightest 
neutralino is (1) bino dominated, (2) higgsino dominated, and (3) v c dominated. For 
each of these cases we consider both the normal and the inverted hierarchical pattern of 
neutrino masses. We show that in these different cases, the ratio of branching ratios of 
certain decays of the lightest neutralino correlates with the neutrino mixing angles. In some 
cases the correlation is with the atmospheric and the reactor angle and in other cases the 
ratio of the branching ratios correlates with the solar mixing angle and in some cases there 
is no correlations at all. Let us now study these possibilities case by case. The interesting 
difference between this study and similar studies with bilinear R-parity violating scenario 



1 25] in the MSSM is the presence of a gauge singlet neutrino dominated lightest neutralino. 
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We will see later that in this case the results can be very different from the bino or higgsino 
dominated lightest neutralino. The lightest neutralino decays in neutrino mass models with 
spontaneous R-parity violation have been studied in ref.Q. 

6.1.1 Bino dominated lightest neutralino 

We will assume that the gaugino masses are unified at the grand unified theory (GUT) 
scale. At the EW scale the ratio of the U (1) and SU (2) gaugino masses are M\jMi = 1:2. 
If in addition, M\ < fx and the value of k is large (so that the effective gauge singlet neutrino 
mass 2kv c is large), the lightest neutralino is essentially bino dominated and it is the LSP. 
First we consider the case when the composition of the lightest neutralino is such that, the 
bino-component I-/V71I 2 > 0.92 and neutrino masses follow the normal hierarchical pattern. 
We have observed that for the bino dominated case, the lightest neutralino (x?) couplings 
to t±-W* pair (where f = e,/i or t) depend on the quantities hi along with a factor 
which is independent of various lepton generations. Naturally, we would expect that the 
ratios of various decay branching ratios such as BR(x° — ► e + W), BR(x7 — ► M + W), and 
BR(x? — ► r + W) show nice correlations with the quantities bf/b 2 with i,j being e, fi or r. 
This feature is evident from fig]?]. Here we have scanned the parameter space of the three 
neutrino Yukawa couplings with random values for a particular choice of the couplings 
A, k and the associated soft SUSY breaking trilinear parameters, as well as other MSSM 
parameters. The trilinear soft parameters A v corresponding to Y v s also vary randomly in 
a certain range. In addition we have imposed the condition that the lightest neutralino 
(which is the LSP) is bino dominated and neutrino mass pattern is normal hierarchical. 



We have checked that the correlations between the ratios of the lightest neutralino 
decay branching ratios and b 2 /b? is more prominent with increasing bino component of 
the lightest neutralino. Note that when (bi/bj) 2 — > 1 the ratios of branching ratios shown 
in fig.0 also tend to 1. We have seen earlier that the neutrino mixing angles #23 and 
6>i3 also show nice correlation with the ratios b 2 Jb 2 and b 2 /b 2 , respectively (see FigsJH 
and H|). Hence we would expect that the ratios of the branching ratios br(^OtVK) anc ^ 
BR{x 7 >e W) show correlations with tan 2 603 and tan 2 #13. These cor- 

relations are shown in fig.||. We have seen earlier (see eq. fl5.12|) ) that with the normal 
hierarchical pattern of the neutri no masses, the atmospheric mass scale is determined by 
the quantity Qf, = ^Jb 2 + b 2 l + b 2 .. Naturally one would expect that the atmospheric and 
the reactor angles are correlated with the £ + W final states of the lightest neutralino 
decays and no correlation is expected for the solar angle. This is what we have observed 
numerically. Here we have considered the regions of the parameter space where the neu- 
trino mass-squared differences and mixing angles are within the 3a allowed range shown 
in eqs.(L_l)-(L3). fig.|| also shows the model prediction for the ratios of branching ratios 
where the neutrino experimental data are satisfied. For our sample choice of parameters 
in fig.||, one would expect that the ratio ^ j^j should be in the range 0.45 to 1.25. 

Similarly, the other ratio BR(x 7 —^e w) . expected in this case to be less 
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Figure 7: Ratio ^ j^j versus ^ plot for a bino like lightest neutralino (the LSP) with bino 

component, | 1 2 > 0.92, where i, j, k — e, /i, r. Neutrino mass pattern is taken to be normal hi- 
erarchical. Choice of parameters are M\ = 110 GeV, A = 0.13, k ~ 0.65, = 300 GeV and to^ = 
400 GeV. Mass of the LSP is 106.9 GeV. The value of the fi parameter comes out to be -228.9 
GeV. 



than 0.07. We can also see from fig.|| that the ratio of branching ratios in the (n + W) and 
(r + W) channels becomes almost equal for the maximal value of the atmospheric mixing 
angle (#23 = 45°). On the other hand, we do not observe any correlation with the solar 
mixing angle #12 since it is a complicated function of of and bf (see eq. (|5.31| )). 

In the case of inverted hierarchical mass pattern of the light neutrinos, the 
coupling is still controlled by the quantities bf. Hence the ratios of the branching ratios 
discussed earlier, show nice correlations with bf/bj (see fig.|9]). However, in this case the 
solar mixing angle shows some correlation with the ratio BR(x 7 —^e w) 

This is shown in fig. 10. The correlation is not very sharp and some dispersion occurs 
due to the fact that the two heavier neutrino masses controlling the atmospheric mass 
scale and solar mass-squared difference are not completely determined by the quantities 
bf and there is some contribution of the quantities af, particularly for the second heavy 
neutrino mass eigenstate. The correlation of the ratio ^|^,o ^ j^j with tan 2 823 shows 
a different behaviour compared to what we have seen in the case of normal hierarchical 
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Figure 8: Ratio * r( *^ Z\ versus tan 2 23 (left), Br( x °-» e wq - h 

& Br ix" >T W) \ J, ^/Br(x° 7 ni W) 2 +Br( X ° >t W) 2 

tan 2 (right) plot for a bino dominated lightest neutralino (the LSP) with bino component, 
I-/V71I 2 > 0.92. Neutrino mass pattern is normal hierarchical. Choice of parameters are same as 
that of fig. 7. 

scenario. This is because in the case of inverted hierarchical mass pattern of the neutrinos, 
tan 2 823 decreases with increasing & 2 /6 2 . One can observe from Figs.|| and 1C that if the 



experimental value of the ratio BR (x 7 > e w ) is <C 1 then that indicates 

a normal hierarchical neutrino mass pattern for a bino-dominated lightest neutralino LSP 
whereas a higher value (~ 1) of this ratio measured in experiments might indicate that 
the neutrino mass pattern is inverted hierarchical. Similarly a measurement of the ratio 
br(^ ir w] can a ^ so ^ ve an indication regarding the particular hierarchy of the neutrino 
mass pattern in the case of a bino dominated LSP. 

6.1.2 Higgsino dominated lightest neutralino 

When one considers higher values of the U(l) gaugino mass Ml, i.e. Mi > /i and large 
value of n (so that the effective gauge singlet neutrino mass 2kv c is large), the lightest 
neutralino is essentially higgsino dominated and it is the LSP. Naturally one needs to 
consider a small value of the coupling A so that the effective /i parameter (ju = 3A^ C ) is 
smaller. In order to look at the lightest neutralino decay branching ratios in this case, we 
consider a situation where the higgsino component in x® is \Nj^\ 2 + | -ZV74 1 2 > 0.90. As in 
the case of a bino dominated LSP, the generation dependence of the X7~^i~W couplings 
comes through the quantities bf. However, because of the large value of the r Yukawa 
coupling, the higgsino-r mixing is larger and as a result the partial decay width of \j in to 



(W + r) is larger than into (W + fi) and (W + e). This feature is shown in fig.ll, where 



the ratios of branching ratios are plotted against the quantities bj /6 2 . The domination of 
BR(x7 — > t + W) over the other two is clearly evident. Nevertheless, all the three ratios of 
branching ratios show sharp correlations with the corresponding bf/b?. In this figure the 
normal hierarchical pattern of the neutrino masses has been considered. As in the case of a 
bino LSP, here also the ratios ppl^o ^ ^1 and BR(x 7 —^e w) show nice 

BR(x" >r W) y/BR(-x° >[i W) 2 +BR(x° 7 >r W) 



correlations with neutrino mixing angles #23 and #13, respectively. This is shown in fig. 12 
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Figure 9: Ratio g^^o ^- j^-j versus ^j- plot for a bino like lightest neutralino (the LSP) with bino 

component | JVy i | 2 > 0.95, where k = e,[M,T. Neutrino mass pattern is inverted hierarchical. 
Choice of parameters are Mi = 105 GeV, A = 0.15, k = 0.65, rape = 300 GeV and m L = 445 GeV. 
Mass of the LSP is 103.3 GeV. The value of the [i parameter comes out to be -263.7 GeV. 



However, in this case the predictions for these two ratios are very different from the bino 



LSP case. The expected value of the ratio rp ,~q_ 



w) ■ 



is approximately between 0.05 and 



AT— ""WO 

0.10 in a region where one can accommodate the experimental neutrino data. Similarly, 
the predicted value of the ratio " ' "' ' 1 



hand, there is no such correlations with the solar mixing angle 612 



is < 0.006. On the other 



Similar correlations of the ratios of branching ratios with bf /6| are also obtained for 
a higgsino dominated LSP in the case where the neutrino mass pattern is inverted hier- 
archical. Once again it shows that the x? decays to (r + W) channel is dominant over 
the channels (e + W) and (// + W) for any values of bf/b? because of the larger r Yukawa 
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Figure 10: Ratio br(x 7 >e w) with tan 2 Qyy (left) plot for a bino dominated 

y/BR(x°—,f, W ) 2 +BR(x° 7 -^T W) 2 y 1 F 

lightest neutralino (LSP) with bino component | | 2 > 0.95. In the right figure the ratio 
gj^u ^ jyj versus tan 2 #23 is plotted. Neutrino mass pattern is assumed to be inverted hier- 
archical. Choice of parameters are same as that of fig. 9. 

coupling. On the other hand, the correlations with the neutrino mixing angles show a 
behaviour similar to that of a bino LSP with inverted neutrino mass hierarchy though with 
much smaller values for the ratios ^1 and BR(x 7 ~* e w ) These 

BR( X %-+T W) ^ B R{x%->fJ. W) 2 +Br(x°->r W) 2 



are shown in fig. 13. Note that the correlations in this case are not very sharp, especially 
with tan 2 0\2- Thus we see that small values of these ratios (for both normal and inverted 
hierarchy) are characteristic features of a higgsino dominated LSP in this model. 

6.1.3 v c dominated lightest neutralino 

Because of our choice of parameters i.e., a generation independent coupling n of the gauge 
singlet neutrinos and a common VEV u c , the three neutralino mass eigenstates which 
are predominantly gauge singlet neutrinos are essentially mass degenerate. There is a very 
small mass splitting due to mixing. However, unlike the case of a bino or higgsino dominated 
lightest neutralino, these u° dominated lightest neutralino states cannot be considered as 
the LSP. This is because in this case the lightest scalar (which is predominantly a gauge 
singlet sneutrino v c ) is the lightest supersymmetric particle. This is very interesting since 
usually one does not get a v c as an LSP in a model where the gauge singlet neutrino 
superfield has a large Majorana mass term in the superpotential. However, in this case the 
effective Majorana mass term is at the EW scale and there is also a contribution from the 
trilinear scalar coupling A k k which keeps the mass of the singlet scalar sneutrino smaller. 
It is also very interesting to study the decay patterns of the lightest neutralino in this case 
since here one can probe the gauge singlet neutrino mass scales at the colliders. 



Before discussing the decay patterns of the lightest neutralino which is v c dominated, 
let us say a few words regarding their production at the LHC. The direct production of v c 



- 32 - 



O 

u 
CO 



0.13r 



0.09 



A 



& 0.05I 



NORMAL HIERARCHY 



LSP HIGGSINO DOMINATED 



0.9 



1.1 



b 2 /b 2 



1.3 



O 0.006 r 



0.004 



'u 0.002 

A 
I 
I 

o 

a O.L 



NORMAL HIERARCHY 



LSP HIGGSINO DOMINATED 

0.02 O04 O06 0.08 

b 2 /b 2 



0.09 r 



NORMAL HIERARCHY 



0.07 



u 

a 



0.05 



A 0.03 



o.<y 



.01 



LSP HIGGSINO DOMINATED 

0.03 O05 O07 0.09 

b 2 /b 2 

e |i 



Figure 11: Ratio jfj^p ^ Yv) versus % pl°t f° r a higgsino like LSP with higgsino component 
(|^73| 2 + l-^74| 2 ) > 0.95, where i,j,k = e,fi,r. Neutrino mass pattern is normal hierarchical. 
Choice of parameters are M x = 325 GeV, A = 0.06, k = 0.65, = 300 GeV and m £ = 400 GeV. 
Mass of the LSP is 98.6 GeV. The value of the fj, parameter comes out to be -105.9 GeV. 



(by v c we mean the v c dominated lightest neutralino in this section) is negligible because 
of the very small mixing with the MSSM neutralinos. Nevertheless, they can be produced 
at the end of the cascade decay chains of the squarks and gluinos at the LHC. For example, 
if the next-to-next-to-lightest SUSY particle (NNLSP) is higgsino dominated (this is the 
state above the three almost degenerate lightest neutralinos) and it has a non-negligible 
mixing with v c (remember that the higgsino-z/ c mixing occurs mainly because of the term 
\v c H\H2 in the superpotential), then the branching ratio of the decay H — > Z + v c can be 
larger than the branching ratios in the 0.W and vZ channels. This way one can produce 
v c dominated lightest neutralino. Similarly, a higgsino dominated lighter chargino can also 
produce gauge singlet neutrinos. Another way of producing v c is through the decay of an 
NNLSP fi, such as f\ — > r + v c . A detailed discussion of these issues is beyond the scope 



of the present paper and we hope to come back to this in a future publication |27|| . 



When one considers higher value of the gaugino mass, i.e. M\ > \i and a small value 
of the coupling n (so that the effective Majorana mass of u c is small, i.e. 2kv c < fi), 
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the lightest neutralino is essentially v c dominated. As we have mentioned earlier, in this 
case the LSP is the scalar partner of v c , i.e. u c . However, the decay of z/ c into v + v c is 
suppressed compared to the decays v c — > + W and v c — > vi + Z that we have considered 
so far. Because of this in this section we will neglect the decay z/ c — > v + v c while discussing 
the correlation of the lightest neutralino (x?) decays with the neutrino mixing angles. 

In this case the coupling of the lightest neutralino (x?) with tj—W pair depends on 
the v c content of x?- Note that the v c has a very small mixing with the MSSM neutralino 
states. However, in some cases the v c dominated lightest neutralino can have a non- 
negligible higgsino component. In such cases the coupling yfy—li—W depends mainly on 
the quantities On the other hand, if X7 is very highly dominated by u c , then the 
coupling X7- £j— W has a nice correlation with the quantities aj. So in order to study the 
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decay correlations of the v c dominated lightest neutralino, we consider two cases (i) v° 
component is > 0.99, and (ii) v c component is > 0.97 with some non-negligible higgsino 
admixture. 



The correlations of the decay branching ratio gR |-o ^ are shown in fig. 14 for the 



cases (i) and (ii) mentioned above. As we have explained already, this figure demonstrates 
that in case (i) the ratio of the branching ratio is dependent on the quantity a 2 /a 2 whereas 
in case (ii) this ratio is correlated with & 2 /6 2 though there is some suppression due to large 
r Yukawa coupling. 

Similar calculations were performed also for the other ratios discussed earlier. For 
example, in fig. 15 we have shown the variations of the ratio ^ ^ as functions of 



2 q2 

-f and rf for the cases (i) and (ii), respectively. The variation with -f is not sharp and 

b 2 

dispersive in nature whereas the variation with rf is very sharp and shows that in this case 
the relevant couplings are proportional to b e ancl b^, respectively. 

On the other hand, in case (i) only tan 2 #23 shows a nice correlation with the ratio 
^ j^j (see fig.lC) and tan 2 6*12 or tan 2 #13 does not show any correlation with the 



other ratio. The non-linear behaviour of the ratios of branching ratios in case(i) is due 
to the fact that the parameters Y u s (which control the aj) appear both in the neutralino 
and chargino mass matrices. The charged lepton Yukawa couplings also play a role in 
determining the ratios. One can also see that the prediction for this ratio of branching 
ratio for case (i), as shown in fig.[l(], is in the range 0.5-3.5, which is larger compared to 
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the bino dominated or higgsino dominated cases for both normal and inverted hierarchical 
pattern of neutrino masses. Also, the nature of this variation is similar to what we see 
with the inverted hierarchical pattern of neutrino masses in the bino or higgsino dominated 
cases. 

In case (ii) none of the neutrino mixing angles show very good correlations with the 
ratios of branching ratios that we have been discussing. However, one can still observe 
some kind of a correlation of the ratio -B-R(x 7 ~^ e w ) with tan 2 6?i o • The 

prediction for this ratio from the neutrino data is on the smaller side (~ 0.07). 
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With the inverted hierarchical neutrino mass pattern, in case (i) one observes a sharp 

BR(x° 7 ^e W) 



correlation of the ratio BR ^ 7 



and«? 



♦e W) 



BR(x° 7 - 



with < 

+T W) a% 



(see fig. 17). The other two ratios BR ^o_ 



^ w) 



.''R( x u >t 1F) ^° no ^ snow ver y snar P correlations with ^f- and |§-, respectively and 
we do not plot them here. However, in case (ii) all the three ratios show nice correlations 
with the corresponding bf/b?. We have shown this in fig. 17 only for b 2 jb\. In this case the 



variations of the ratios of branching ratios with neutrino mixing angles are shown in fig. 18 



For the case (i) , only tan 2 #13 shows certain correlation with the ratio of branching 
ratio shown in fig.|l8| (right), but we do not show it here. 

Finally, we would like to mention that in all these different cases discussed above, the 
lightest neutralino can have a finite decay length which can produce displaced vertices in 
the vertex detectors. Depending on the composition of the lightest neutralino, one can have 
different decay lengths. For example, a bino-dominated lightest neutralino can produce a 
displaced vertex ~ a few mm. Similarly, for a higgsino dominated lightest neutralino, 
decay vertices of the order of a few cms can be observed. On the other hand, if the lightest 
neutralino is v° dominated, then the decay lengths can be of the order of a few meters. 
These are very unique predictions of this model which can, in principle, be tested at the 
LHC. 

7. Summary and Conclusion 

In this work we have studied a supersymmetric model in detail where the observed pat- 
tern of neutrino mass-squared differences and mixing angles are obtained with the help of 
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three standard model gauge-singlet neutrino superfields, which simultaneously solve the 
H problem of MSSM. The additional terms in the superpotential and the scalar potential 
include R-parity violating interactions involving these gauge-singlet neutrino superfields. 
The vacuum expectation values of the singlet sneutrinos give rise to effective Majorana 
mass terms for the singlet neutrinos, as well as a /i-term, both at the electroweak scale. 
This model was introduced in ref.|| and some phenomenology was discussed for a single 
gauge-singlet neutrino superfield. The spectrum and parameter space of this model, with 
three gauge singlet neutrino superfields, were discussed in ||. We have performed a de- 
tailed and extensive analysis of this model in the neutrino and neutralino sector with the 
inclusion of three generation of gauge singlet neutrino superfields along with the associ- 
ated interaction terms. The neutrino mass matrix is obtained because of the electroweak 
scale seesaw mechanism involving the gauge-singlet neutrinos and the mixing between the 
MSSM neutralinos and the neutrinos. We have done a thorough and systematic study of 
the neutrino mass matrix both analytically and numerically and tried to identify the rele- 
vant parameters which crucially control the bilarge pattern of neutrino mixing angles. We 
show that even with a flavour diagonal structure of the neutrino Yukawa coupling matrix, 
two large and one small mixing angles can be generated in this case. Both the normal and 
inverted hierarchical pattern of neutrino masses can be obtained with different hierarchies 
of the neutrino Yukawa couplings. Because of the presence of the neutrino- neutralino mix- 
ing, it is in general difficult to obtain a degenerate mass spectrum of the neutrinos and we 
do not consider this possibility in this paper. 

We have also looked at the scalar sector of this model and wrote down the neutral 
scalar, pseudoscalar and charged scalar mass-squared matrices of this model. The allowed 
regions in the model parameters are obtained which satisfy certain constraints in the scalar 
sector. For example, absence of tachyons in the scalar squared-mass eigenvalues puts severe 
constraints in the relevant parameter space. With these choices of parameters, satisfying 



- 38 - 



the scalar sector constraints, we have tried to fit the global three flavour neutrino data with 
both normal and inverted hierarchical mass spectrum of the neutrinos. Since this model 
involves several free parameters, we have made a few simplifying assumptions which makes 
this model more predictive. 

Perhaps the most interesting part of this whole analysis is the study of the phenomenol- 
ogy of the lightest neutralino which can be the LSP in some cases or it can be the NLSP. 
The decay patterns of the lightest neutralino may provide additional information to find 
out more about the neutrino mass patterns and mixing angles. We have considered three 
different scenarios where the lightest neutralino can be either a bino-dominated one or a 
higgsino-dominated one or it can be mostly a gauge-singlet neutrino with very little mix- 
ing with the other states. We have also discussed briefly the production mechanism of the 
gauge-singlet neutrino dominated lightest neutralino at the LHC. The study of the decay 
pattern and the production mechanism of the lightest neutralino in this last mentioned 
scenario is extremely important because it will help in probing the mass scale and the 
properties of the gauge-singlet neutrinos at the LHC. The presence of the gauge-singlet 
neutrino dominated lightest neutralino can also distinguish this model from the usual bi- 
linear R-parity violating model of generating neutrino masses and mixing. An important 
test of this model, as a super symmetric solution to the observed neutrino mass patterns 
and mixing, can be performed by measuring the ratios of the decay branching ratios of the 
lightest neutralino in the final states involving a charged lepton and a W-boson. We have 
shown explicitly that these ratios of branching ratios have certain correlations with the 
neutrino mixing angles which depend on the nature of the LSP as well as on the pattern of 
the neutrino mass hierarchies considered. The study of the higgsino dominated LSP case is 
also very important because it can provide information about the /x parameter which is de- 
termined in terms of the vacuum expectation values of the gauge-singlet sneutrinos. Thus 
one may have information about the gauge-singlet neutrino mass scale and its coupling, 
from the decay pattern of the higgsino dominated LSP. Collider phenomenology of this 
model at the LHC is very rich both in the fermionic and the scalar sector. For example, 
pair produced lightest neutralino at the LHC give rise to the final states fj,fj,WW, ttWW 
and firWW with a certain ratio for their production rates nicely correlated with tan 2 #23- 
These production rates also depend on the dominant component of the lightest neutralino 
as well as on the different hierarchical patterns of the three light neutrino masses. Another 
important testable prediction of this model is the measurement of displaced vertices in the 
decay of the lightest neutralino. This decay length can vary in the range of a few mm to 
~ 1 meter depending on the nature of the lightest neutralino. We hope to come back to 
these issues in a future publication [27]. 
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A. Scalar mass squared matrices 

In this appendix we present the details of various scalar mass squared matrices. For 
the convenience of the reader let us repeat that the scalar mass squared matrices are now 
8x8, considering all three generations of sneutrinos (both doublet and singlet) and charged 
sleptons. These enhancements are essentially due to the mixing of neutral Higgs bosons 
with both the doublet and singlet sneutrinos and the mixing of charged Higgs with the 
charged sleptons. 

A.l Neutral scalar mass squared matrices 

The decomposition of various neutral scalar fields in real(7£) and imaginary(Z) parts are 
as follows 

H i = H m + iHiii 
"I = Hn + ^fcx> 

(A.l) 



Only the real components get VEVs as indicated in eq, Q2.3| ). 

The entries of the scalar and pseudoscalar mass-squared matrices are defined as 

i T,/r2\aB _ I ^ & ^neutrals 

{Mp) - { 2^4 } > 

(M 2 )a p = (l^^L), (A.2) 

where 

4>j = Hii, H^j-, v%x-> Vk2, 

= H m ,H2Ti,i>kTz,i>k-R.- ( A -3) 

Note that the Greek indices a, (3 are used to refer various scalar and pseudoscalar Higgs and 
both doublet and singlet sneutrinos, that is H®, H®, v\\-,Vk-, whereas k is used as a subscript 
to specify various flavours of doublet and singlet sneutrinos i.e., k = e, fj,, r in the flavour 
(weak interaction) basis. 

A. 1.1 CP-odd neutral mass squared matrix 

The basis for CP-odd or pseudoscalar mass-squared matrix is 

*l = (H% z ,H% E ,i%r,i>ra). (A.4) 
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The pseudoscalar mass term in the Lagrangian is of the form 

r mass , , = $^2^ 

'"pseudoscalar ^pj»ip^_K, 



(A.5) 



where A4p is an 8x8 symmetric matrix. Using eq.( [2.7D , (2.8), and eq.( |2.9[ ), the independent 
elements are given by 



(M 2 P ) H ^ 



1 



j \ ik / j j i 



(M|) 


n 2T n 2T 


(Mj.) 


n \X v mI 


{Ml) 




(Ml) 




{Ml) 


ttQ r, 
2X u rnl 


{Ml 




{Ml 




{Ml 





i,3,k 
1 

V'2 



j \ l,k J i,j i 



7, C 



2^<V - E(^ F ^) im ^ + (^AA) m ^i, 

i,j,k j 

_ 2 J2 K jnm C j + 4$ju™iu$ + p m p n + h nm v\ + (m? c ) nm - 2^2{A K K) inm ^, 
j j i 

2 s £u^Y^v 2 - Y™Yj>i + r > m + l*viY? m ~ A m r>i - (A u Y v ) nm V2, 
j » 

(A.6) 



where 



X n x m + ^ywyjm 



(A.7) 



We have checked that one eigenvalue of this 8x8 matrix is zero corresponding to the 
neutral Goldstone boson. 

A. 1.2 CP-even neutral mass squared matrix 

The basis for the CP-even or scalar mass-squared matrix is 



^5 - ( H lKi H 2Ki Kill Vnn)- 
The scalar mass term in the Lagrangian is of the form 



(Ai 
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(A.9) 



where M 2 S is an 8 x 8 symmetric matrix. The independent entries using eq.(p.7|), (|2.g|), 
and eq.(|2.9|) are given by 



1 



E A ^2 ( E^' fc ^ ) + E AVw 2 2 + + E^am^ 

j \ ifc / j j i 



(M 



J)^2-R.^2T. 



1 

^2 



-E^ E^V** ) -E^w'^+E^a)^ 



h3 



(Mj)^V^ = -2^A J <% 2 + 2 f iv 1 X m - A m E r >i " ^ m ~ (^AA) m w 2 , 



j 

(M 2 s ) H ™* c m n = 2E<V + 2A m /i« 2 + 2^y« m r> 2 + EOW™^ " (^aA)™^, 

j » i 

(Mj)^V^ = 2E*TV^2 + YY^^ k vtu% - 2 lg u m v 2 + EOW™^, 

(M s ) mn = 22^k j C +42^ c J U C J +p p +h v 2 + (m c) + 22_^{A K n) u i: 

j j i 



YjWv\ + r™r n c + lg C v dnm + 2 lg u n v m + (m|) 



2 Nmn 



(A.10) 



where eq.( |A.7|) has been used. 

A. 2 Charged scalar mass squared matrix 

The charged scalar mass squared matrix considering all three generations of both left 
handed and right handed charged sleptons is also an 8 x 8 matrix. For the sake of complete- 
ness, here we give the expressions for various elements of the charged scalar mass-squared 
matrix. The elements of the charged scalar mass-squared matrix are defined as 



{ c) ~ { 2 dd ca dd cP 



where 



n l i n 2 ' e Kk> e Lk 



(A.ll) 
(A.12) 
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The basis for charged scalar mass-squared matrix is 

*f = (H+,H+,e+ n ,et n ), 
and the charged scalar mass term in the Lagrangian is of the form 

pmass _ jfy— T k/1 

charged scalar ~ JVl C^C> 



(A.13) 



(A.14) 



where M.q is an 8 x 8 symmetric matrix. The independent elements of M? c using eqs,(|2.7|), 



2.8j ), and eq. (|2.9]) are given by 



(Mq) HiHi 



{Mq) HiH2 



(ii4) HiS «" 



(M 2 C ) H ^ 



(M^) H2iR 



(Mcf Rn ~ eLr - 



(M, 



1 



i,j,k i 



1 

«2 



3 3 i 

i 

r, 2 

i 92 /V"^ 2 , 2\ 



u l 1 



7, C, 



I -J 



i i 
3 » 

n 2 



'•J 



-\xY e v 2 + (A e Y e ) vi, 



Ln e Lm — n 

' c ' c 



ffl; 



si. 



+ ^y e ^V e "^f + la t v 5 mn - ^ v 5 mn + fv m v n + (ml) 



2 \mn 



(A.15) 



As mentioned earlier, we have computed the eigenvalues of the charged scalar mass- 
squared matrix numerically and ensured that seven of its eigenvalues are positive and there 
is a charged Goldstone boson. 

B. Feynman rules 

In this appendix we will study the relevant Feynman rules for the LSP decay calculations 
28. 29]. The required Feynman rules are 

1. Neutralino-neutralino - Z°, 



2. Chargino-neutralino - 



- 43 - 



B.l Neutralino-neutralino-Z and chargino-chargino-Z ,7 

For neutralinos the following relations between mass and weak eigenstates are very useful 



P L B° 




PrB° 


= PrNhxI 


p l wI 


= P L N* 2 xl 


PrWI 


= P R Ni2Xi, 


PlH 3 


= PlN* j+2 xI 


PrHj 


= P R Nij + 2Xi, where 


P L v k L -- 


= PlN* k+7 xl 


Pr4~- 


= PRN it k+7Xi, 


Pl4 


= P L N* k+4 xl 



Prv c R = P R N iM4 xt where k = 1,2,3, (B.l) 
with i varies from 1 to 10 and 




(B.2) 



In terms of the four component spinors Xi f° r charginos, the following relations between 
mass and weak eigenstates are very useful 



PlW = P L V£xi, 
P R W = P R UaXi, 
P L H = P L V* 2 Xi, 
P R H = P R U i2 Xi, 
P L V = P L V* j+2 Xi, 

P R P = P R U itj+2 Xi, (B.3) 

where j = 1,2,3, and i varies from 1 to 5. 

So in terms of physical or mass eigenstates of charginos and neutralinos the required 
interaction terms are as follows 



z°xx 



r 



9 



+ 



cos 9w 
9 



2 cos 9w 



z^xh* 1 



O'^Pl + O'^Pr 



x 



Xj: 



(B.4) 
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where 



Og = -ViiVfr - -V l2 V* 2 + 5ij sin 2 9w, 

O'ij = -UtxUji - ^U* 2 U j2 - ^U* k+2 U jM2 , +5ij sin 2 9 W , 



o" R 



1 N i3 N* 3 + \N^N* A - }rN iM7 N* k+7 , 



2 w JO 2 
-Of*,fc = l J 2,3. 



In deriving eq. flB.q ) unitary properties of U and V matrices has been used. 

,n 

(ig/2 C 9w ) ^ [o£ (1 - 7 5 ) + 0* (1 + 7 5 )] 





Y 11 

-V; 



(i ff /2 cor [c^(i- 7 5 )+of (1+7 5 )] 



(B.5) 





it 



ig/2)^ [o£(l- 7 5 ) + 0*(l+7 5 )] 



where Cg w = cos 9w 

Figure 19: Feynman rules for interaction of neutral and charged gauge boson with charginos and 
neutralinos. 

B.2 Chargino-neutralino - 

Now in terms of physical chargino and neutralino states the relevant interaction term is 



'W-x+x 



+C 



w+ x -x°> 



where C w +^-^o is the hermitian conjugate of C w -^+^o and 

C w - t +x o = gW^-f [Ot 3 P L + 0§P R ] x} 



(B.6) 



(B.7) 
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with 



0% = N i2 V* x - ^N l4 V* 2 , 

0§ = N* 2 U n + -^=N* 3 U j2 + l=N* k+7 U jtk+2 , 



where k = 1, 2, 3. 

The Feynman rules are shown in fig.@ The matrices Ogs Og 1 , C^ L , Of, R and 0£, 



Ofj are defined by eq. (|B.5| ) and eq.( |B.§| ), respectively. 
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